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ABSTRACT 

Analyses  show  that  the  central  place  in  the  mathematical  description  of  the  spatial-temporal  dynamics  of  one  or 
several  interacting  populations  take  the  task  of  research  of  interaction  of  the  migration  population  processes  with 
demographic.  Population  models  that  take  into  account  only  the  demographic  processes  are  well  known  and  fairly  well- 
developed  [3,4].  This  so-called  point  model,  the  main  assumption  which  is  the  assumption  of  "infinitely  fast"  stirring  of 
individuals  in  a  given  area.  This  assumption  is  true  if  modeled  area  is  quite  small  compared  with  the  mean  free  path  of 
individuals,  or,  equivalently,  with  a  radius  of  individual  activity  [3].  If  this  provision  violated,  in  population  models  must 
account  migration. 

KEYWORDS:  Biological  Population,  Nonlinear  Model,  Parabolic  Type  Equation,  Self,  Similar  Solutions,  Lower  and 
Upper  Solutions 

INTRODUCTION 

The  simplest  and  most  widely  used  in  the  present  situation  is  the  hypothesis  of  randomness  "walk"  of  individuals 
in  space.  This  assumption  allows  substantiating  of  using  as  a  tool  of  modeling  equations  of  reaction-diffusion  type,  where 
as  a  reaction  part  was  used  the  right  part  of  the  point  models,  and  diffusion  coefficients  (mobility  of  individuals)  are 
assumed  to  be  constant. 

In  the  framework  of  these  models  it  is  possible  to  explain  such  effects  as  wave  propagation  of  population  at 
settling  the  area  and  existence  of  complex  spatio-temporal  dynamics  of  population. 

Consider  in  QT  =[0,T]x  RN  generalized  reaction  -  diffusion  task  of  Kolmogorov  -  Fisher  type  in  the 
following  form 
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which  describes  the  process  of  biological  populations  with  double  nonlinearity,  which  diffusion  coefficient  is 
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Equation  (1)  is  a  generalization  of  the  simplest  diffusion  model  for  the  logistic  model  of  population  growth  [1].  It 
can  be  regarded  as  the  equation  of  nonlinear  filtration,  thermal  conductivity  in  the  presence  of  simultaneous  influence  of  a 
source  and  acquisitions. 

This  equation  in  the  case  of  m=0,  k(t)=  k  >0  constant  qas  suggested  (as  well  as  the  cubic  equation  instead 
quadratic  nonlinearity  in  the  right  side)  Fisher  [3]  as  a  stochastic  model  of  the  deterministic  version  of  the  gene  in  a 
favorable  distribution  of  diploid  populations.  He  examined  in  detail  the  equation  and  obtained  some  useful  results. 
Heuristic  and  based  on  genetics  conclusion  of  the  equation  have  also  led  by  A.  N.  Kolmogorov,  I.  G.  Petrovskii  and 
N.S.Piskunov,  classical  work  [4]  which  served  as  the  basis  for  a  more  rigorous  analytical  approach  to  the  Fisher  equation. 

Properties  of  the  solution  of  the  Cauchy  problem  and  the  solution  of  equation  (1)  in  the  case  of  a  homogeneous 
environment,  when  (1)  m  —  0  and  k(t)=k  -  constant  studied  in  detail  by  many  other  authors  (see  for  example  [3],  where 
you  can  find  links  to  other  work).  Suppose  that  in  ( 1)  m  —  0,  k(t)  =  k  -constant.  In  this  case,  for  speed  of  propagation  of 
the  wave  with  wave  solutions  of  (1) 


authors  of  [4]  were  rated  an  estimate  C  >  . 

With  regard  to  the  properties  of  solutions  of  the  initial  task  for  the  equation  (1)  is  unknown  and  so  it  is  interesting 
to  trace  the  evolution  of  the  process  of  reaction  diffusion  in  heterogeneous  environment  and  explore  the  impact  of 
heterogeneity  and  case  of  dependences  of  the  reaction  coefficient  from  the  time,  i.e.  k=  k(t). 

This  paper  is  devoted  to  the  study  of  properties  of  solutions  of  the  task  (1),  (2).  Clarified  conditions  on  parameter 


estimation  of  a  solution  of  the  task  (1). 

We  construct  self-similar  equations  for  (l)-(2)-more  simple  for  research  equations.  We  construct  self-similar 
equation  by  the  method  of  nonlinear  splitting  [1] 

Below  we  propose  an  algorithm  based  on  the  splitting  of  the  original  equation  of  parabolic  type  with  which 
assesses  estimations  of  lower  and  upper  solutions  of  the  task  (l)-(2). 


Df"-cf'  +  f(l-fp)  =  0 
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where  take  the  place  finite  speed  of  movement  of  a  wave.  Proved  also,  bilateral 


I.  ALGORITHM  OF  SPLITTING  OF  THE  PARABOLIC  TYPE  EQUATION 


We  solve  the  equation 
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(3) 


and  then  search  solutions  of  the  equation  (1)  in  the  form 


u(t,  x)=u  (t)w(r,  x) 


(4) 


Then  supplying  (4)  into  (1)  we  have  again  equation  of  the  form  (1) 
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where k,(t)  =  k(t)u'n(p-2)+m  =  k(t)um<"-l) . 

only  difference  is  that  instead  of  the  t  stands  r(t)  =  J 
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Hence,  if   mp  —  m  =  0  equation  (5)  has  exactly  the  form  (1),  so  again,  we  get  an  equation  of  the  form  (1) 
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however,  instead  of  a  variable  by  the  time  t  stand  the  function  v(t)  =  j[u(t)]    f  dt. 
Now  consider  equation 
dw 


dr 


=  VQ  Vw  ~  \p'2  Vw  ) 


(7) 


Equation  (7),  which  is  called  the  etalon  for  (1),  has  six  types  of  self-similar  solutions,  one  of  which  has  the  form: 
W(T,x)  =  f(a     Z=\X\/TUP  (8) 


where  /(£)  satisfies  the  equation 
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at  ^  <  (a  I b)(p  1)1  p  .  From  (6)  it  is  clear  that  if  selected  as  the  comparable  function  z(t,  x)  =  u(t)w(r,  x)  and 
W  <  1  in  QT  ,  where  W  is  solution  of  (6),  u(t)  given  by  (3),  it  will  be  a  lower  solution  of  the  original  task,  if 
U0  (x)  >  z(0,  x)  .  This  result  can  be  improved  by  applying  again  algorithm  of  nonlinear  splitting  applied  above. 

Now  after  putting  w(t,  x)  =  f(^),    £  =1  X  I  /  TUp ,  equation  (6)  will  turn  into  a  self-similar  equation 
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Now  if  we  construct  a  function  z(t,x)  =  u(t)w(T,  x)  ,  where  w(t,  x)  -  upper  solution  of  the  etalon  equation 

(7),  for  which  runs  conditions  Az  <  0  in  QT  and  U0  (x)  <  z(0,  x)  ,  then  function  z(t,  x)  will  be  an  upper  solutions  of 
the  task  (1),  (2).  Here 


Aripov  M.  M.  &  Muhamediyeva  D.  K. 


du  d 

Au  =  

dt  dt 


D 


V 


du' 


dx 


p-2  \ 

du"'  ] 


dx 


+  k(t)u(l-up). 


On  these  considerations  are  based  algorithm  for  the  construction  of  upper  and  lower  solutions  of  the  Cauchy  task 

(1),  (2),  the  global  solvability  of  the  initial  task  and  definition  of  the  critical  value  of  the  parameter  @  ,  where  solution  of 
the  task  becomes  unlimited,  and  changes  view  of  the  asymptotic  representation. 

,    ,  2  =2,  D={(t,x):  t>0, lxk/(0}     ,     i(t)>0  f>n 

Let's  introduce  notations  ^     ^  lv     ;  WJ   ,  where  'w^u  at  1  - u  -  continuous 


function. 


Takes  place  the  following  lemma  of  comparison  of  solutions  [3,5]. 

Lemma.  Let  U^'X^  —  ^  -  generalized  solution  of  the  task  (1),  (2)  and  functions  u±(^x^  sucn  that  U±  ~  ^ 
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and    in    D    runs    an  inequality 
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Functions  U+ ,  U_  are  called,  respectively,  the  upper  and  lower  solution  of  (1),  (2). 
By  direct  calculation  we  can  ensure  that  the  function 
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where  Cj  <  0,    O  =  a      ,  is  generalized  solution  of  the  equation 
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where  constant  a  >  0  . 
H.  ESTIMATION  OF  SOLUTION  OF  THE  TASK 

Theorem  1.  Let  in  (l)-(2)      «0(jc)  <  u+(0,x),  x  e  7?^     =  x/r"/'  and 


(12) 
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Then  in  Q  there  is  a  global  solution  of  (1),  (2),  for  which  we  have  the  estimate  u(t,  x)  <  W+  (f,  x)  ,  where 
u+(t,x)  =  u(t)f0(g),  (13) 


here  w(0>  /o     )  given  by  formulas  (3)  and  (11). 
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Suggested  that,  /?  >  1 ,  0  <  U0(x)<  1 ,  0 <  k(t)e  C(0,oo)  .  For  solution  of  this  task  is  fair 

Theorem  2.  Let  0  <  w0  (x)  <  1   X  e  R  N  .  Then  for  the  solution  u(t,  x)  of  the  task  (1),  (2)  in  Q  takes  place  an 


estimate 


u(t)(TQ  +tyNe  V4D,)  <u(t,x)<er>  (T0+ty"e 


jk(t)dt 
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where  u  (V)  -solution  of  the  equation  (1)  without  diffusion  part: 
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and  constant  7"0  >  1 . 


Study  various  properties  of  the  solution  of  the  Kolmogorov-Petrovsky-Piskunov-Fisher  task  with  diffusion  with 
reaction  applied  to  the  study  of  biological  populations  attracts  the  attention  of  many  researchers,  as  it  will  be  of  interest 
from  the  point  of  view  of  mathematics  as  a  nonlinear,  which  has  an  important  application  task.  Therefore  it's  important 
generalizate  this  task  on  different  cases. 

In  particular,  it  has  not  been  studied  in  case  of  a  heterogeneous  environment  (diffusion  coefficient  is  a  function  of 
the  spatial  variable),  reaction  coefficient  depends  on  time  or  has  a  more  complicated  character.  Consider  the  task  of 
reaction  diffusion  of  Kolmogorov  -  Fisher  type 


du  _  d 
dt  dt 
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(17) 


describing  the  process  of  diffusion-reaction  in  heterogeneous  (anisotropic)  environment  with  heterogeneity 
coefficient  x"  ,  where  O^neS. 

Propose  an  algorithm  for  constructing  upper  solutions  of  the  Cauchy  problem  (16), (17)  by  the  method  of 
nonlinear  splitting  and  etalon  equations. 

For  the  solutions  of  the  Cauchy  task  (16),(17)  fair  following  theorem: 

Theorem  3.  Suppose  that  in  (16)  performed  following  conditions 

U(0,X)  <  Z(0,X)  ,  where  z(t,x)  =  U(t)  •  /(£)  , 


U  (f )  is  a  solution  of  the  equation  Ut  —  kit)u  (1  —  U  ^  )  , 
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Then  the  problem  (16),(17)  is  globally  solvable  and  for  the  solution  in  the  field  of  Q  take  place  following 

estimation  u{t,x)  <  z(t,x) . 

I  \n 

Presence  of  coefficient  X    in  (16)  makes  difficult  of  study  the  task  (16),  (17)  in  terms  of  building  self-similar  or 

some  partial  solutions,  as  well  as  numerical  modeling.  Turning  to  the  new  variables,  assuming  ll(t,X)  =  w(t,(p(x)) 

and  choosing  (p(.x)  ,  reduce  (16)  to  a  form  which  in  the  case  j3  >  1  previously  more  detail  studied  by  many  authors  [1], 
PL  [5]. 

HI.  NUMERICAL  SCHEMES  AND  METHODS  OF  SOLVING 

To  solve  numerically  the  task  (1)  -  (2)  in  the  let's  build  uniform  grid  03 h  by  x  with  step  h: 

G)h  =  {x{  -  ih,    h>0,    i  -  0,1,...,  n,    hn  =  b], 

and  temporary  grid 

WT  =  \tj  =  jr,    t  >  0,    j  -  0,1,...,  m,    r  m-  t\,  T>0. 

Replace  the  task  (18)-(20)  by  using  balance  method  with  implicit  difference  scheme  and  receive  differential  task 
with  error  0(h2  +  t)  . 

Now  we  present  the  results  of  numerical  experiments.  Iterative  process  is  based  on  the  method  of  Picard  and 

Newton. 

Results  of  computational  experiments  show  that  all  iterative  methods  are  suitable  for  the  constructed  schemes. 
Below  listed  results  of  numerical  experiments  at  k(t,  x)  —  k(t)  in  a  heterogeneous  environment  in  two-dimensional  case, 
for  different  values  of  parameters  included  in  the  equation  (Fig.l): 

t  g  [0,1];  xt  e  [-6,6];  x2  e  [-6,6];  nXj  =  30  ;  nXi  =  30  ;  nt  =  100 
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Values  of  Parameters 


Results  of  Computational  Experiment  at 
the  Initial  Moment  of  the  Time 


Results  of  Computational  Experiment 
at  the  Final  Moment  of  the  Time 


m(p-\)  +  n-  \  >0 
w  =  0.7,  P  =  \.5 

n=3.1,p=4.1 
eps  =  10  ~3 

Average  It  1=1, 43 
Average  It2=  1 
Average  It=2,43 


time(l) 


time  (100) 


m(p-\)  +  n-  \  =0 

m  =  0.5,  yff  =  1.5 

n=l,p=0.5 

eps  =  10  3 

Average  It  1=2, 19 
Average  It2=l,51 
Average  It=3,70 


time(l) 


time(lOO) 


-0.2 


.0.1 


m(p-l)  +  n-l<0 
m  =  0.1,  £  =  1.5 

n=l,p=0.1 

eps  =  10  ~3 
Average  It  1=2,02 
Average  It2=l,21 
Average  It=3,23 


time(l) 


time(lOO) 


Figure  1:  Results  of  Numerical  Experiments  in  the  Two-Dimensional  Case 
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